Spectral properties of finite laser-driven lattices of ultracold Rydberg atoms 
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We investigate the spectral properties of a finite laser-driven lattice of ultracold Rydberg atoms 
exploiting the dipole blockade effect in the frozen Rydberg gas regime. Uniform one-dimensional 
lattices as well as lattices with variable spacings are considered. In the case of a weak laser coupling, 
we find a multitude of many-body Rydberg states with well-defined excitation properties which are 
adiabatically accessible starting from the ground state. A comprehensive analysis of the degeneracies 
of the spectrum as well as of the single and pair excitations numbers of the eigenstates is performed. 
In the strong laser regime, analytical solutions for the pseudo-fermionic eigenmodes are derived. 
Perturbative energy corrections for this approximative approach are provided. 

PACS numbers: 32.80.Rm, 32.80.Ee, 37.10.Jk 



I. INTRODUCTION 

In the past decades, experimental investigations of ul- 
tracold atoms have reached an unprecedented control 
over the motional degree of freedom as well as the in- 
teraction of the atoms. One of the intriguing systems 
in this field are ultracold Rydberg atoms since they ex- 
hibit strong long-range interactions of isotropic or even 
anisotropic character. In fact, the interaction-induced 
level shift of collective states featuring multiple Rydberg 
excitations can be large enough to exceed the excitation 
laser linewidth thereby preventing the excitation of fur- 
ther Rydberg atoms. This dipole blockade effect was pre- 
dicted by theory one decade ago along with a proposal 
how it could be exploited for the realization of fast two- 
qubit quantum gates (TJ [2] . Since then, a wide range of 
more general quantum information processing applica- 
tions has emerged, putting Rydberg atoms in the focus 
of many theoretical as well as experimental efforts (see [3] 
and references therein) . Experimentalists have found the 
dipole blockade effect [?HB] in gases of alkali atoms and 
recently for just two individually trapped atoms [51 110] . 
The coherence of the observed effects was demonstrated 
[lH - TH] and the resulting collective Rabi-frequency [T5] 
has been observed. 

The difficulties often encountered in working with a 
large number of atoms are related to the spatial struc- 
ture of their gaseous samples. Because of the varying 
particle density typically encountered in a trap and the 
spatially varying intensity profile of the excitation laser, 
no unique collective Rabi-frequency is encountered. Sev- 
eral theory groups worked to circumvent these limitations 
by performing simulations of the excitation dynamics in 
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unordered samples looking for signatures of the coherent 
dynamics |15H17| and predicting the formation of crys- 
talline structures |18j . In [TS] it was proposed to use a 
laser detuned from the single atom resonance to actively 
produce Rydberg atom pairs with the specific distance 
that corresponds to an interaction induced energy shift 
equal to the detuning. This effect was observed indirectly 
by measuring the interaction-induced ionization rates as 
a function of the interaction time |2U] . In spite of these 
successes, it remains clear that a spatially ordered Ryd- 
berg gas is highly desirable. To this end, both theorists 
as well as experimentalists have concentrated on iden- 
tifying a stable trapping mechanism for Rydberg states 
that allows a strong confinement in magnetic |2 11425] as 
well as optical [26 -28 traps. While the creation of lattice 
traps is straightforward nowadays in the optical regime, 
ongoing experimental effort is put into creating arrays of 
magnetic traps |29l 150] . 

First theoretical investigations of the Rydberg excita- 
tion in structured ultracold atomic gases [3T] consider 
a ring-shaped lattice of ground state atoms, concentrat- 
ing on the case where the laser coupling to the Rydberg 
state is weak in comparison to the next-neighbour Ryd- 
berg interaction. In the opposite regime of a dominant 
laser coupling, the same authors demonstrated that the 
system permits the formation of fcrmionic collective ex- 
citations |52"] . 

In the present work, we consider a finite one- 
dimensional lattice of ground state atoms that are co- 
herently excited to the Rydberg state via a two-photon 
laser transition. We provide a thorough investigation of 
the spectral properties of this system for single-spaced 
as well as variably-spaced lattices. Different parameter 
regimes are covered. The model Hamiltonian that forms 
the basis of our investigations is derived in Sec. [IT] In Sec- 
tion |III A"] the weak laser regime is analyzed where we find 
that a multitude of system states with well-defined ex- 
citation properties are adiabatically accessible from the 
ground state. Among these are the crystalline states that 
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have been previously discussed in [15]. In Section IIIB 



we focus on the strong laser regime where we achieve an 
approximate description through the XY-model of the 
antiferromagnetic spin chain |33) . We provide analyti- 
cal solutions for the pseudo-fermionic eigenmodes as well 
as the perturbative energy corrections that are necessary 
due to the approximate nature of the description. The 
intermediate regime between a weak and strong laser 
coupling is briefly treated numerically in Sec. |III C| In 
Sec. |IV A| we consider the case of multiple lattice spac- 
ings and describe in Sec.[lVB]a specific spatial setup that 
allows the emulation of the crystalline state transitions of 
a larger sized lattice. The appendix contains a derivation 
of the perturbative energy corrections in the strong laser 
regime. 



II. THE MODEL 

Starting from a quite general model including all in- 
dividual atoms contributions, one can derive a simpli- 
fied effective Hamiltonian that describes the system in 
terms of a lattice of two-level systems. To this end, 
the internal structure of an individual atom is described 
in terms of a ground state g, an excited (Rydberg) 
state e, and an intermediate state m that we include 
in order to allow for an experimentally realistic excita- 
tion scheme via a two-laser setup: g < — > m < — > e. 
After transforming into a rotating frame of reference 
and a subsequent coarse-graining in time, an effective, 
static Hamiltonian is obtained and the intermediate m 
state can be removed from the state space by an adi- 
abatic elimination if the excitation lasers are strongly 
detuned with respect to this level |341 155] , After this, 
the two-laser setup is modelled by an effective g < — > e 
coupling and the description at each lattice site k can 
be further reduced by employing the superatom states 

^| ff f\ ff f\...,ef\...,^) and 

, g^l/i respectively Here Nk is the 
number of atoms on site k. This means that either all 
atoms of a given site are in the ground state or they 
symmetrically share a single excitation [12] . This re- 
duction is justified only if the confinement to a single 
lattice site is sufficiently strong such that all atoms are 
located within a distance of each other that is smaller 
than the dipole blockade radius and also much smaller 
than the lattice constant. The symmetrization of the ex- 
cited superatom state presupposes that both the effective 
coupling strength of the two-photon transition and its de- 
tuning from resonance are identical for all atoms within a 
given site. By considering solely the internal level struc- 
ture of each atom, the frozen Rydberg gas regime is pre- 
sumed where the atomic centre of mass motion can be 
neglected on the timescale of the coherent excitation dy- 
namics. 

The above considerations lead to our model Hamilto- 
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Here, the operators <J^ k \i € {x, y,z\ act on the su- 
peratom located at site k and take on the usual Pauli- 
matrix form when expressed in the local superatom basis 
SW ■= {\e) k ,\g) k }- The excitation number operators 
may also be expressed in terms of the Pauli-operators, 

The system is therefore 
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formally equivalent to a spin- 1/2 lattice with interactions 
Vki- In this picture, the contributions due to the laser 
are similar to the interaction of the spins with an exter- 
nal magnetic field with a (local) component aligned 
with the spins and a perpendicular component fi'*' . 

Hamiltonian ([l]) contains three different contributions. 
First, the laser coupling of each single atom's ground 
state to the excited state is given by 
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where fiw := ^/N^o denotes the collective Rabi fre- 
quency for the superatom state at site k [T2]. Because of 
the number of atoms Nk contributing to the excitation 
dynamics, fi^ is enhanced by a factor of y/Nk compared 
to the single atom Rabi frequency £Iq. Second, the part 
which describes the site-dependent laser detuning, 
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corresponds to an energy gap (in the effective RWA pic- 
ture) of A 1 -^ between the excited state and the ground 
state of the superatom at site fc. Finally, the Rydberg 
interactions read as pairwise interactions between each 
two sites k, I with an interaction strength that depends 
on the spatial separation of the sites, 
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Each summand is non-zero only if both affected sites are 
in the excited state. We only consider repulsive inter- 
actions, i.e., Vfcf > 0, which are common for interacting 
Rydberg atoms in their ns-state for a wide range of prin- 
cipal quantum numbers n |36| . 



The contributions (24) may be classified by two sep- 
arate criteria: By locality, i.e., whether or not they act 
non-trivially on more than a single site, and by their 
simultaneous diagonalizability. The laser contributions, 
i.e., the laser coupling Hl and the laser detuning Ho 
are local, while the Rydberg interactions _ff int are by 
definition non-local. On the other hand, the Rydberg 
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interactions and the laser detuning operator commute. 
[Hjj, Hi n t] = 0, since they can both be expressed in terms 

(k) 

of a z -operators and the identity. We thus identify two 
interesting limiting parameter regimes: In the weak laser 
coupling regime the Hamiltonian is dominated by its di- 
agonal contributions: the Rydberg interactions and the 
laser detuning. Alternatively, in the strong laser regime 
the Hamiltonian is dominated by local operators. As 
demonstrated in [35] for a ring lattice, this allows for 
an approximation of the system by an XY-model. Un- 
less stated differently, we will restrict ourselves to global 
laser parameters (flw, A^) — ► (ft, A). 

We describe the system in terms of the 'canonical' 
product basis 

Sn ■= {\sis 2 ■ ■ -s N ) , s k e{e,g}}, (5) 

since these states are directly accessible in experiments. 
Moreover, our Hamiltonian is already diagonal in this ba- 
sis except for the laser coupling part Hl ■ Calculating the 
matrix elements of our Hamiltonian H for the canonical 
product basis Sn yields a sparse matrix. Specifically, 
it is straightforward to show that the number of non- 
zero matrix elements is given by (N + 1)2 = D log 2 2D, 
where D = #5jv = 2 N is the dimension of the state 
space. To simplify our notation, we define the canoni- 
cal product ground state and the fully excited state as 
| G) := | gg . . . g) and \E) :— \ ee . . . e) , respectively. 

III. SINGLE-SPACED LATTICES 

We start by considering lattices with a single lattice 
spacing a. In this case, the interaction potential is given 

by 

We will usually assume a Van-der-Waals interaction po- 
tential, i.e., n = 6. The restriction to global laser param- 
eters yields the final Hamiltonian 

k k l>k 1 1 

This Hamiltonian is invariant under reflections at the 
centre of the lattice. As was done for the ring lattice 
in [32J we designate the corresponding unitary operator 
for this symmetry as 1Z and define it via its action on the 
local Pauli operators, 

for ii = x, y, z. The product ground state is invariant 
under reflections 1Z\G) = \G) and since the full basis 
can be constructed by means of the ground state and the 
Pauli-operators, this completely determines the form of 
1Z. Clearly, reflecting the system twice should leave it 



unchanged and hence TZ^ = TZ^ 1 = 1Z. The eigenvalues 
of 1Z are thus given by ±1. 

For a weak laser coupling |ft| <C V\, |A|, the diagonal 
contributions to the Hamiltonian dominate. In this case 
the laser coupling leads to a small off-diagonal perturba- 
tion. Alternatively, in the case of weak Rydberg interac- 
tions the Hamiltonian can be mapped approximately to 
an XY-model Hamiltonian. In the following, both lim- 
iting regimes are discussed in detail and analytical for- 
mulas describing the excitation spectra are derived. For 
the intermediate regime, where neither of the above con- 
ditions is fulfilled, numerical simulations reveal the full 
spectrum. 

A. Weak Laser Regime 

We start by considering the weak laser regime, i.e., 
assuming |ft| < Vi,|A|. The Hamiltonian m is conve- 
niently divided into two parts, H = H + W, grouping 
together the laser detuning with the next-neighbour Ry- 
dberg interactions to give the dominant contribution 

H = 1 -Aj2*i k) ' + ViY:ni k H k+1 \ (8) 

k k=l 

while the perturbation consists of the laser coupling as 
well as the long range Rydberg interactions, 

N-x N-d 
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which we have rewritten as a sum of contributions for 
a given separation d. Introducing the operators for the 
total excitation number, 

N 1 - N 

^. = E»? ) = 5E^ ) + T' (10) 
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and the next-neighbour excitation pair number, 

JV-l 

N ee =yni k H k+1 \ (11) 

k=l 

we can rewrite the unperturbed Hamiltonian as Hq = 
A(N e ~N/2) + V x N ee . For a given state \S) = 
\s-\_S2 ■ ■ ■ sn) € Sn from the canonical product basis 
[cf. ([5j], the unperturbed energy eigenvalue is thus given 

by 

E(S) = A[N e (S)-N/2] + V 1 N ee (S), (12) 

where N e (S) and N ee (S) denote the eigenvalues of the 
operators N e and N ee for the state | S) . They can read- 
ily be obtained by counting the number of excitations 
. . . e . . . and excitation bonds . . . ee . . . present in the se- 
quence S — sis 2 ■ ■ ■ sn- There are several important ob- 
servations to be made at this point: 
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1. Ha is linear in N e and N ee . Since N e and N ee are 
diagonal, this linear energy relation holds for their 
eigenvalues as well. 

2. In general, there are multiple canonical product 
states of equal (N e ,N ee ) which are consequently 
always degenerate with respect to the unperturbed 
Hamiltonian. 

3. Depending on the specific ratio of A and Vi, the 



simple form of ( 12 ) already suggests that additional 



degeneracies between states of different (N e ,N ee ) 
are possible. Since degeneracies in the RWA picture 
correspond to resonant laser couplings in the non- 
rotating frame, our result simply states that we can 
tune the laser to resonantly excite multi-particle 
states, as one would intuitively expect. 

We define Dj^(N e ,N ee ) to be the dimension of 
each (N e , -/V ee )-subspace "HtN t ,N ee )- The calculation of 
Dn(_N b , N ee ) is possible through a combinatorial analy- 
sis. First, note that for a given state \S) = \ sis 2 ■ ■ ■ sjv) 
the number of next-neighbour excitation pairs N ee (S) 
is fully determined by the total number of excitations 
N e (S) and the number of excited domains ee . . . e within 
\S) which we denote by d e (S). Whenever d e (S) is equal 



to one. N P „ takes on its maximal value N Pl 



N„-l. 



'ee\ de = 1 

For each additional domain N ee decreases by one if N e 
remains fixed. Hence, the following relation holds for any 



d e (S) =N e (S)-N ee (S). 



(13) 



We can now calculate Dn by analyzing the number of 
possibilities of how to construct appropriate sequences 
's\S2 ■ ■ • sjv'. For N e — we can only have IS") = 
\gg,..g) — \G). hence we assume N e ,d e > 1. We must 
distribute N e excitations across d e < N e domains which 
leads to a factor of \ j"Zi) = (n -jv^-i)- Now, for any 
such division of the excited states into domains, we must 
count the number of ways how to distribute these do- 
mains of excited atoms across the lattice such that there 
is always at least one ground level site between two ex- 
cited domains. The number of ground level sites is given 
by N g = N — N e . The domains of excited sites can thus 
be inserted at N g — 1 positions between two 'g '-characters 
or at the two positions at either end of the lattice. Hence, 
there are N g + 1 positions across which we distribute d e 
excited domains and we must multiply the above result 
by a factor of ( N ° H +1 ) = ( T tr Nc N +1 ). Together, we find 



D N (N e > l,iV ee ) 
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and D N {0,N ee ) = 5 Neefi - 

The spectrum of Hamiltonian (|8|, i.e., the full set of 
unperturbed energy eigenvalues as a function of a vary- 
ing ratio A/Vi exhibits points of high degeneracy, see 
figure [T] These occur only at specific rational values of 
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FIG. 1: Energy spectrum for N = 13 lattice sites and a 
constant laser coupling Q, = 0.05Vi. At specific rational values 
of A/Vi many energy levels are almost degenerate. Depending 
on A/Vi the energetic ground state in the RWA frame is given 
by (1) the canonical product ground state, (2) the alternating 
state | ege . . . ge) , and (3) the fully excited state. 



A/Vi. In principle, these can be calculated from (12 1 
if one takes into consideration the combinatorially possi- 
ble combinations of N e and N ee . For A = the lasers 
are tuned to atomic resonance. In this case the state 
| gg , . . g) is degenerate (within the effective RWA picture) 
with all states that lack neighbouring excitations which 
corresponds to dipole-blocleqed states. However, for non- 
zero detuning the state \gg . . . g) can also be brought 
to degeneracy with other states. Physically, this cor- 
responds to a situation in which the lasers resonantly 
couple | gg . . . g) to states containing pairs of neighbour- 
ing excitations. Note, however, that due to the local 
nature of the laser coupling Hamiltonian Hl, these cou- 
plings require the presence of intermediate, in general 
off-resonant states. A state from a given (N e ,N ee ) sub- 
space can be coupled via Hl to states from (N e ±1, N ee ), 
{N e ± 1, N ee ± 1) and (N e ± 1, N ee ± 2). 

A particular example is provided by A/Vi = 
— 1/2: In this case \gg...g) is degenerate with 
states containing exactly two neighbouring excitations 
{I e-eg ...g), \geeg . . . g) , . . . \g . . . gee)}. The general con- 
dition for a degeneracy between two canonical product 
states |5i) and j^) is given by 



A iV ee (5i)-iV ee (S 2 ) 
Vi N B (Si) - N e (S 2 ) ' 



(15) 



For iV e (5i) ^ N e (S 2 ) and N ee (Si) ^ N ee (S 2 ) the equa- 
tion can only be solved if the ratio A/Vi is rational be- 
cause the eigenvalues of N e and N ee are integral. A com- 
binatorial analysis reveals that the quantum numbers N e 
and N ee obey the following constraints in addition to be- 
ing integral: 

N e = => N ee = 0, (16) 
1 < N e < \N/2] ^0<N ee <N e -l, (17) 
\N/2] <N e =>2N e -N-l<N ee <N e -l. (18) 
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These render it difficult to provide a concise, general 
quantitative analysis of the possible degeneracies. 

However, to gain a better understanding of the spec- 
trum and its degeneracies we can visualize the combina- 
torial possibilities for (N e , N ee ) and try to recognize gen- 
eral principles that are independent of the lattice size. To 
this end, figure [2] illustrates all possible pairs (N e ,N ee ) 
for N — 13. For a fixed ratio A/Vi, (12 1 can also be 



interpreted as describing the lines of constant energy in 
the (N e ,N ee ) plane. These lines can be constructed by 
drawing a line of slope A = — A/Vi through each valid 
(N e ,N ee ) point. A degeneracy occurs whenever two or 
more points lie on the same line. The number £ of re- 
maining non-degenerate energy levels depends sensitively 
on the ratio A/Vi. The minimal such number £ mm , i.e., 
the highest degeneracy is attained for A/Vi = — 1. This 
is true independently of N since it follows from the shape 
of the distribution of possible (N e ,N ee ) pairs. How- 
ever, the actual value of £ m ; n depends on the system 
size and is given by £ = [AT/2] + 1. The maximal value 
of £, on the other hand, is determined by the possible 
(N e ,N ee ) pairs for a fixed number of lattice sites and 
reads £ max = 1 + \N/2] 2 . Results of a systematic cal- 
culation of the number of non-degenerate energy levels 
for many different ratios A/Vi are shown in Figure [3] 
for N — 13. When moving away from the point of the 
highest level of degeneracy (A/Vi = —1) the number 
of non-degenerate energy levels increases from its mini- 
mal value £ m i n = 8 and rapidly reaches its maximal value 
£max = 50, except for the points of integer A/Vi for which 
pronounced local minima of £ are encountered. 

In the following, we proceed by analyzing the most 
relevant parameter regimes in more detail. We are espe- 
cially interested in the level structure of the energy band 
that includes the canonical product ground state | G) or 
the fully excited state \E). For a strongly blue detuned 
laser A/Vi ~ —2, for example, we find that the fully ex- 
cited state \E) is degenerate with all other states that 
fulfil N ee = 2N e - N - 1 > 0; the latter are given by 
the right hand side of the red 'triangle' of states compat- 
ible with the constraints, cf. figure [2] These states mini- 
mize the next-neighbour Rydberg interaction energy for 
a given number of excitations N e . Another interesting 
parameter setting is given by A/Vi ~ —1/2. This leads 
to the resonant creation of isolated pairs of neighbouring 
excitations. Specifically, the subspace that includes the 
product ground state | G) is characterized by the relation 
2N RK = N R . 



1. The Full Blockade Regime 

The full blockade regime is attained by restricting both 
the laser detuning and the laser coupling to be very small 
compared to the next-neighbour Rydberg interaction en- 
ergy. The relevant subspace that includes the canonical 
product ground state is then approximately (neglecting 
the mixing due to a non-zero f2) given by the 0-eigenspace 



of N ee , i.e., all states that contain no neighbouring exci- 
tations as given by the lowest row in figure [2] Note that 
this implies N e < [jV/2] . 

For non-zero detuning, the full-blockade subspace is 
energetically separated from the remaining state space if 
the detuning fulfils |A| < Vi/N. Introducing the gener- 
alized excitation pair number operators 
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we may rewrite our full Hamiltonian as 



A - O N w " 1 1 



(20) 
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For a positive detuning A > 0, the canonical product 
ground state is clearly also the ground state of the diag- 
onal part of Hamiltonian (20 1. However, at A = 0, even 



when considering the part of the Rydberg interactions 
beyond nearest neighbours, this state becomes degener- 
ate with all (N e — 1, N ee — 0) states. For small negative, 
i.e., blue detuning, these states form a degenerate ground 
state energy level, up until A/Vi = V(A — 1)" at which 
point they cross with the N e = 2 state | eg . . . ge) . Further- 
decreasing the detuning leads to a succession of ground 
states with an increasing number of excitations 2 < N e < 
\N/2~\ . In general the ground state energy level is degen- 
erate, except when (N— l)/(N e — 1) — i m j n > 2 is integral: 
In this case, it is possible to evenly distribute the exci- 
tations across the lattice, hence, there exists a unique 
minimal energy crystal state \ [N e ]) := \eg . . . geg . . . ge) 
where the spacing between neighbouring excitations is 
given by Z m i n - This state's eigenvalue for the generalized 
excitation pair operators is given by 

Am i fat i\ _ /(^ - I M for « = e N < 

eeU el) -\Q otherwise. 

(21) 

The corresponding diagonal matrix elements of the 
Hamiltonian accordingly read 



([Ne]\H\[N t 



- N/2) + A £ "' 



in x — , „„ 

mm q —i y 



Vi 



(22) 



A(N e -N/2) + j-^(N e -l) (23) 



A(N e - N/2) + 



\n+l 



(n — iy 



(24) 



When (N — l)/(N e — 1) is non-integral, the minimal en- 
ergy configurations are realized by states with excitations 
separated by at least two different spacings. This case is 
more complex to analyze in full generality. 
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FIG. 2: Lines of constant energy (yellow) in the (N e , N ee ) plane for N = 13 and for different slopes (a) A/Vi = —1/2 and (b) 
A/Vi = — 1, respectively. The line thickness indicates the dimension of the associated energy eigenspace. The dashed red lines 
illustrate the constraints [see (16|17|18| >]. The size of an X-marker indicates the dimension D]si(N e ,N ee ) of the corresponding 
subspace. The minimal number of non-degenerate energy levels is attained for A = — Vi, cf. subfigure (b). Note that in this 
case the canonical product ground state (N e = 0, N ee — 0) is energetically isolated. 
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FIG. 3: The number £ of non-degenerate energy levels for 
different ratios 0^ for N = 13. For any N, the minimum 
number is attained for A = — Vi. In our particular case we 
find £ m i n = 8 and £ max = 50, respectively. 



If both N e , N e - 1 e N> 2 are divisors of N - 1, the 
crystal states | [N e + 1]) , | [N e ]) feature an avoided level- 
crossing for a well-defined detuning A^ e+1 : 



Vi 



N — 1 J I ^ k n 

7 Kk=l 



(25) 
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In the thermodynamic limit N — > oc 
/ € (0,1/2] of the exact expression 
suits of Weimer et al. are reproduced, 
f n [(n + l)C(n) + 0(1/N e )} [37!. For our modest lattice 
sizes this limit is not properly realized, but we can see 
that the scaling of the transition detunings with the ex- 



NJN 
j25| the 
-A/Vi 



citation fraction / = N e /N is determined by the scal- 
ing of the interaction potential with the distance. For 
van der Waals interactions (i.e., n = 6) this implies that 
the transition detunings for small excitation fractions are 
very small. 

Whenever one is interested in the crystal states | [N e ]), 
it is favourable to choose a lattice size N for which many 
different | [N e ]) —> | [N e + 1]) transitions can be realized. 
For example, one can verify that N = 13 allows for the 
transitions | [2]) ->■ |[3]) |[4]) ->■ |[5]). The next can- 
didate of an ideal lattice size is given by N = 61 which 
allows for the additional transitions | [5]) — > | [6]) — > | [7]). 



In Section IV we present an alternative approach how to 
realize these transitions by means of a non-uniform lat- 
tice with a specific pattern of lattice spacings. Figure [4] 
illustrates the succession of crystal states for the above- 
mentioned case of = 13. An interesting aspect about 
this phenomenon is that although the locations of the 
level crossings, i.e., the specific values of the detuning at 
which they occur, are determined by the scaling of the 
interaction potential with the distance, the existence of 
such crystal ground states already follows from any kind 
of long-ranged repulsive interactions. 



2. Fully Excited Lattices 

We conclude our discussion of the weak laser regime 
by considering a strongly blue detuned laser A/Vi ~ —2. 
In this case, we find that the fully excited state \E) is 
degenerate with all other states that fulfil N ee — 2N e — 
N — 1 > 0. In Figure [2] these states make up the right 
hand edge of the triangle of states compatible with the 
constraints. 

As we have seen in figure [l] for a strongly blue detuned 
laser, A < — 2Vi, the fully excited state \ ee. . .e) actually 
becomes energetically favoured in the RWA frame. Be- 
cause of the long-ranged Rydberg-Rydberg interactions 
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FIG. 4: (a) Detail of the energy spectrum near the ground 
state for N = 13 lattice sites in the full blockade regime for 
small detunings and very small laser coupling Q,/Vi = 0.005. 
The red li ne i ndicates the ground state. At specific detunings, 
given by ( 25 1 , avoided crossings are encountered where the 
ground state changes from one crystal state to another, (b) 
Site-resolved local excitation probability of the ground state 
as a function of the same laser detuning as in subfigure (a). 
We see from right to left a sequence of states featuring an 
increasing number of evenly spaced excitations. White regions 
correspond to unit Rydberg excitation probability, whereas 
dark regions denote the ground state. 



one finds a succession of states with an increasing num- 
ber of Rydberg excitations for increasing negative detun- 
ing. This trend continues until the completely excited 
lattice is encountered. Hence, the A/Vi « —2 transition 
concludes what was started at the A = degeneracy: a 
transition from the canonical product ground state to the 
maximally excited state with all sites occupied by Ryd- 
berg excitations. The A/Vi ~ —2 transition is illustrated 
in figure [5] 



3. Resonant Pair Creation 

As mentioned before, another interesting degeneracy is 
realized for A/Vi = —1/2. From the degeneracy condi- 



FIG. 5: Same as in figure [4] but for detunings close to 
A/Vi = — 2. One encounters the transition from the alter- 
nating excited state | ege . . . ge) (which is the final state is 
figure [4| to the fully excited state \E). 



tion ( pL5[ ) one finds that in this case the state | G) (N e — 
N ee = 0) is degenerate with all other canonical prod- 
uct states whose quantum numbers satisfy N e = 2N ee . 
Figure [6^ a) presents the spectrum at the point of degen- 
eracy for N = 13. The maximal number of excitations 
compatible with both the combinatorial constraints and 
N e = 2N ee is given by the maximal integer satisfying 
N e < |(iV - 1). For N = 13 this is N e = 8 and conse- 
quently N ee — 4. Correspondingly, in subfigures gb) and 
(c) the projections Vo,o and "Pg,4 are shown for the two 
distinct states that evolve from and into the canonical 
ground state \G), respectively. Vn c N ee is the projector 
onto the subspacc 'Hn,,,n cc containing N e Rydberg exci- 
tations and N ee Rydberg pairs. As expected, one finds 
a transition of the adiabatic eigenstates from the canon- 
ical ground state to the pair excited states when going 
through the region of degeneracy. At the crossing itself, 
further subspaces contribute to the investigated states as 
7^0,0 + ^8.4 7^ 1- The main such contribution stems from 
the 7^6,3 subspace while the tIa.2 an d %2,i subspaces only 
contribute marginally. 
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FIG. 6: (a) Detail of the energy spectrum including the 
canonical ground state \G) for detunings close to A/Vi = 
-1/2, N = 13, and for Q/Vi = 0.1. This specific detuning 
favours the excitation of Rydberg pairs. In subfigures (b) and 
(c) the projections Vo,o and Vsa are shown for the two dis- 
tinct states marked by the dash-dotted green and dashed red 
line. 



B. Strong Laser Regime 

As for the weak laser regime, when the laser energy 
contributions are very large in comparison to the next- 
neighbour Rydberg interactions, the system again be- 
comes analytically treatable. Neglecting all long-range 
interactions, we rewrite the interaction part in this 
regime as 



H, 



fe=l 



Vi 



N 

fc=i 



(*0 



4 



[<x«+aW]+(iV-l)| 



(27) 



The total Hamiltonian can then be separated according 
to H = Ha + H' + Hi, into a local part 



(N) 



an interaction part H' 

boundary term Hf, = — ^j- ( cr)t' + cr^'J. Defining 

the pseudo laser coupling Cl := \j£l 2 + (A + V\) 2 ^> 
Vi, the local part can be re-expressed as Hq := 
|^Sfc=i^^ where we have introduced the notation 



^-^i£{x ,y, 



(k) 



The rii are given by 





(29) 



The pseudo angle 9 G [0, 2ir] satisfies the relation 



cos = — = — , sin = — . 



(30) 



Note that H contains no information about the spatial 
arrangement of the lattice sites. It is symmetric under 
any permutation of the lattice sites. Therefore, it can 
be diagonalized by rotating each site k independently in 

(k) 

spin-space about the <7y axis with the unitary transfor- 
mation 



N 



)exp 



fc=i 



o v 



It follows that 



Ulna^Ue, = cos{8 - 0'W fc > 



sin y - 



9'Vl fc) . 



(31) 



(32) 



Choosing 9' = 6, the transformed local Hamiltonian is 
given by 



H = U%H U e = if>I>« + (TV - 1)^. (33) 



k=l 



From Hamiltonian (33), the spectral features can be 



straightforwardly deduced. Its eigenstates are given by 
products of the rotated single site basis states 



cos - e 
2 1 



'fc+ Sm 2 Iff)*" 



sin 2 l e )fc + cos 2 



(34) 
(35) 



The transformed product ground state and the fully ex- 



cited state read 



G\6 



U e \G) = | ;0), and 

£-;6»^ = U g \E) = | + H h;0), respectively. Due to 

the full permutation symmetry of Hq the associated en- 



ergy depends only on the number N + of 
present in a given product state, 



factors 



1 N 

H Q = -J2{^ ) + (A + V 1 )o^ 



fc=i 



l)f . (28) 



1 w 



k=l 



WgN+Ug, 



(36) 
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FIG. 7: Energies as a function of the pseudoangle 6 for 
Vi = 0.1S1. The deviations from the constant energy level 
behaviour are strongest for f2 = (corresponding to 8 = 
0, 7T, 271"). They are always positive, which is due to the fact 
that the perturbing Rydberg interactions operator is strictly 
positive. 



As pointed out previously, 6' = 9 ss 0, ir corresponds to 
an almost vanishing laser intensity, i.e., we recover the 
weak laser regime (cf. Sec. Ill A I. Consequently, all off- 

(k) 

diagonal terms (i.e., those featuring at least one eri ' op- 
erator as a factor) vanish due to the factor of sin 9. Hence, 
within the vicinity of these values, the energy level struc- 
ture can be fully understood in terms of our previous 
discussion of the weak laser regime. The conceptually 
opposite case is realized for 9' = 9 ss ir/2,3ir/2. In this 
case the laser coupling is maximal, its magnitude being 
given by the pseudo laser coupling which we assume to 
be much larger than the Rydberg interactions and the 
laser detuning, |f2| = O 3> V\ = —A. In the following 
we fix 9 — ir/2; the case 9 — 3tt/2 can be treated equiv- 
alently. As mentioned before, this closely corresponds to 
the system analyzed in [38 . The main differences are 
given by our requirement that A = — Vi and our open 
boundary conditions. Therefore we will present in the 
following discussion only the major differences. 



where N + = | J2k=i [ cos ® a 



sm( 



ai k) + 1]. The dom- 



2 L~ik=\V 

inant part Hq of the Hamiltonian in the strong laser cou- 
pling regime can be expressed in terms of N + and N + 
as 



and 



H = n{N+ - N/2) + (N - 1) 



H = n(N+ - N/2) + (N - 1) 



V 



(37) 



(38) 



For 9 = tt/2 our interaction Hamiltonian H' = 
"4 1 J2k=i fJ a:' C ' cr i' C+1 ' ) assumes a purely off-diagonal form. 
The transformed boundary term accordingly reads Hj, = 

■4 L [cx L ' > + <?4 ]• Clearly, H' is not diagonal in the rotated 
basis, but it is possible to further separate it into a part 
that commutes with Ho and one that does not. This is 



achieved by rewriting = a 



( k) 

■a_ , leading to 



respectively. The eigenvalues of N + are k = 0,1, ... ,N 
and for each eigenvalue < k < N the eigenspace has 
dimension (^) due to the number of possibilities to dis- 
tribute k excitations across a lattice of length N. Hence, 
the eigenvalues of H depend only on Q. Moreover, its 
spectrum is symmetrical with respect to E = for Vx = 
and is given by 



-{N — 2/c)f2/2, 



0,1, 



,N. 



(39) 



This fact is illustrated in figure [7] where we confirm the 
constant spacing between the energy bands. The slight 
variation of the energy levels with 9 stems from the Ry- 
dberg interactions which are incorporated in H' and Hi,. 

To understand the substructure of the energy levels at 
other values for 9, we must focus our attention on the 
interaction part H' . Applying the unitary transform for 
9' = 9 to the interaction Hamiltonian yields 



H' = ^008*9 J2°i k) °i k+1) 

k=l 



fc=l 



Vi . 

— sin tt cos 
4 



N-l 



fc=i 



(40) 



H' = ^J2 [^4 k+1) +a {k) 4 k+1) 



N-l 



k=l 
N-l 



V 

4 



Er (fe) (fe+i) . (fe) (fe+i) 
\(T + a + + <J_ cr_ 



k=l 



H[ + H' 2 . 



(41) 



H[ consists of operators which conserve the total num- 
ber of excitations N + and therefore commutes with 
the unperturbed Hamiltonian, i.e. 



fr (fe) (fc+1) 



fr (k) (k+1) 
Hq,C7_ <T + 



for any k € {1, ... N - 1} . This 



implies that we can find a simultaneous eigenbasis for Ho 
and H[ . As we will show later, the remaining part of the 
perturbation is fully off-diagonal in this basis. Hence, 
for our parameter regime 3> V\ , A the dynamics of 
any initial state that belongs to a certain N + eigenspace 
is influenced by the off-diagonal part H' 2 only at second 
order 0(Vf/fi). Formally, this can be demonstrated by 
the application of quasi-degenerate perturbation theory 
[39 . We thus neglect H' 2 in the following and proceed by 
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diagonalizing the resulting Hamiltonian 

X T N-l 



Defining new fermionic operators r\j :— X^fcLi RjkCk, 
the Hamiltonian decouples into N independent 
'fermionic' modes: 



fe=i 



A' 



Vi 



k=l 



ft^af/' - iVfi/2 + (JV - 1) 



(42) 
(43) 



Vx 
4 



fc=i 



As already pointed out in [3T] , ( 43 1 is the familiar Hamil- 
tonian of the XY-Model [33J , which can readily be diago- 
nalized by the introduction of fermionic ladder operators 



(fc) i 
Cfe := a_ 'e 



(-i) fe -v«n^). (44) 



These non-local operators and their adjoint operators 
obey fermionic anti-commutation rules 



{c k ,ci} = {cl,c}} = 0, {c fc ,cj} = d fc 



fc/- 



(45) 



By repeated action of different c k on the ground state 
G^ we can construct a basis for the state space which, by 

construction, is also an eigenbasis of N + (and therefore 
Hq). It is also straightforward to confirm the following 



relations: 



(fc) (fe) f 
aX'&L> = c k c k , 

(fc) (fc+i) f 
a+(T- '=c k c k+1 , 



-c k c 



t 

fc+1 



c i+i cfc - 



(46) 
(47) 
(48) 



Hence, the XY-Hamiltonian, which is a quadratic form 
in the a± operators, becomes a quadratic form in the 
fermionic operators, 



N N-l 

H xy = 2ClJ2 c{c k + -f { c l Ck + 1 + c l+i Ck 



(49) 



k=l 



NQ/2 + {N - 1) 



fc=l 

Vx 



which, by introducing the real, symmetric matrix M — 
(Mjk)^ k=1 , we may rewrite as 



N N 



4 9 = ^X! c)M jk c k - NQ/2. 

j=l k=l 



(51) 



The elements of M are given by Mj k = £18 j k + ^ (<5j.fc+i + 
fij+i.k)- By diagonalizing M — R T AR, Aj k — Xj5j k with 
an orthogonal matrix R, the Hamiltonian is further sim- 
plified. This corresponds to a principal component anal- 
ysis of the quadratic form in the fermionic operators. 



JV 



H xy = J2 X kvUk ~ NQ/2 + (iV - 1) 



Vi 



(52) 



fc=i 



G 



The system's ground state is still given by 

| ) and the orthogonal transformation of the op- 
erators preserves the anti-commutator relations. There- 
fore, all other eigenstates can now be constructed by the 
successive application of ri\. operators for different k since 
each mode can only be occupied by a single excitation 
[(Vk) 2 ~ 0]- The matrix elements of R and the trans- 
formed operators are given by 

Rjk = 



sin 



n + i v^ + i 

K N 



N ■ 



fc=l 



jkn 

M "' JV + l'"- 



(53) 
(54) 



A straightforward calculation reveals that the eigenvalues 
of M are given by \j = fl + ^cos(-^j-). Hence, our 
resulting Hamiltonian reads 



N 
fc = l 



S I H cos 

2 WV+1 



-NQ/2 + (N-l) 
(55) 



Vi 



The fermionic excitation number operators rj^rjn appear- 
ing in this Hamiltonian have only and 1 as eigenvalues. 
Neglecting momentarily the term — NQ/2+ (N— 1) the 
spectrum is thus given by 



E. 



kik 2 ...k[ 



2n ^ 

n=l 



fc„7r 



le{l,...,N},ki < k 2 <••• < h g{1,2 



,N} 



(56) 



\kik 2 ...h) := vl 1 vl 2 ---vl l 
total excitation number operator obeys 

N N 



G 



(50) The corresponding basis states are given by 

Note that our 

(57) 
(58) 



ivr ( fc ) ( fc ) ' 



k=l 



k=l 



N 



N 



= Y RnkRmkVnVm = X) ^ 



fe,n,m— 1 n— 1 

and counts the 77-mode excitations: 

N + \kik 2 ...ki) =l\kik2.-.ki) 



(59) 
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FIG. 8: Energy spectrum in the strong laser regime, showing 
the ground state, the iV+ = 1 subspace, and partially the 
N+ — 2 subspace for N = 12. Without the effect of the 
perturbation H'2, the energy levels would scale linearly with 
Vi for fixed Q. Here, we see that for large perturbations, the 
linear scaling is slightly violated. For the ground level we have 
also plotted the analytical result of the unperturbed energy 
level (green, dash-dotted line) as well as the result of second 
order perturbation theory (red, dashed line). 



This confirms that the neglected second part of the 
original perturbation H'2 = ^ Y^k=i \. a -i ' 



.C fe ) (J ( fe + 1 )j j g f u Hy off-diagonal in our final eigenbasis, 
as it maps any basis state with I excitations into a super- 
position of states containing I ± 2 excitations. Hence, for 
Cl Vi the off-diagonal perturbation only contributes to 
the spectrum at second order 0(V~i /Cl) and to the eigen- 
states at order 0(Vi/fl), justifying our former approach. 

For the ground state and the first excited states we 
have explicitly calculated the resulting second order en- 
ergy corrections (see appendix). Figure [8] presents the 
numerically calculated spectrum of H as well as the ana- 
lytically calculated results for the ground state. A com- 
parison of the analytical and numerical results for the 
energy eigenvalues of the N + — 1 subspace {ry[. | G) , k — 
1,2,..., N} is provided in table [I] In general, a very good 
agreement between the analytical and numerical results 
is obtained, demonstrating that Hf, and H' 2 only lead to 
minor corrections. 

For completeness, we give a brief discussion of the sym- 
metry properties of the eigenstates. A detailed calcu- 
lation reveals that the transformation properties under 
reflection are 



^v.^ = (-ir +i ^(-i) JV+ = (-i) 



iff. 



(60) 



G), its eigenvalue 



For a single fermionic excitation rf k 
with respect to the reflection symmetry operator is cor- 
respondingly given by (— l) fc+1 . Hence, only states cor- 
responding to an odd k = 1, 3, 5, . . . belong to the same 
symmetry eigenspace as the canonical product ground 

state. For a general state ry ki rj]. 2 ■ ■ ■ ry k G) the eigen- 



TABLE I: Comparison of the N+ — 1 energy levels calculated 
analytically from the unperturbed Hxy Hamiltonian with re- 
sults obtained via second order perturbation theory and the 
exact numerical results for N = 12 and Vi/Cl = 0.1. 



k 


R (0) 


B (o) , AE (2) 


jrinum 


1 
1 




A £81 




9 


4 R81 


4 RRZ, 
_i ±.uou 


4 RS^ 


Q 
O 




4 fiQI 
-*±.uy 1 


4 RQ9 


4 


-4.697 


-4.701 


-4.700 


5 


-4.707 


-4.709 


-4.711 


6 


-4.719 


-4.723 


-4.722 


7 


-4.731 


-4.733 


-4.734 


8 


-4.743 


-4.746 


-4.746 


9 


-4.753 


-4.756 


-4.757 


10 


-4.762 


-4.767 


-4.766 


11 


-4.769 


-4.774 


-4.774 


12 


-4.774 


-4.778 


-4.778 



value is given by (_i)fci+fc2+...fc ( -M(i+l)/2_ 

In concluding this section, we remark that the results 
obtained here only apply in a rotating frame of reference 
due to the RWA picture. The expectation value with 
respect to a given state | ip) of any operator O diago- 
nal in the canonical product basis S can be evaluated as 
usual: (O) = (ip\0\i/j)- Any off-diagonal operator, on 
the other hand, must first be transformed to the RWA 
frame and therefore becomes explicitly time dependent. 
Within the weak laser regime discussed in the previous 
section, most of the relevant observables are diagonal in 
the canonical product basis and are thus not affected by 
the RWA-picture. For the strong laser regime eigenba- 
sis, the situation is more complicated and needs to be 
inspected individually. 



C. Intermediate Regime 

In the intermediate regime, where our parameters 
A, fl, Vi can be of the same magnitude, we cannot sepa- 
rate the energy scales due to Rydberg interactions and of 
the laser contributions. We will therefore limit ourselves 
to the discussion of some numerically obtained results. 
As in the previous section, we fix the relative scale of 
Rydberg interaction energies and laser contributions b 
varying the pseudo angle 9 for fixed and V\. Figure 
provides two examples of the resulting spectrum, namely 
for Vi = ft/2 and Vi = 2ft, respectively. For 9 = 0, tt, 2tt, 
the Hamiltonian is diagonal in the canonical product ba- 
sis iS/v and the spectra can be fully understood from our 
discussion of the weak laser regime. For intermediate 
values of 9, however, our analytical results for the strong 
laser regime break down because the mixing of the en- 
ergy bands due to the neglected part of the interactions 
and the boundary term eventually prevail. Accordingly, 
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FIG. 9: Energy spectrum obtained in the intermediate regime 
as a function of varying 6 while keeping the Rydberg inter- 
action energy and the pseudo laser coupling constant, (a) 
Vi/Cl — 1/2 and (b) Vi/fi = 2. The clear separation of the 
energy bands as seen for the strong laser regime in figure [7] is 
lost. 



in figure [9] the corresponding energy bands cannot be re- 
solved anymore and only the lowest and highest energy 
states remain well separated. 



IV. NON-UNIFORM LATTICES 

In this section we generalize our setup to allow for lat- 
tices with variable lattice spacings. Due to the spatial 
dependence of the inter-site interactions this introduces 
new energy scales into the spectrum. To account for the 
different lattice spacings, we define the spatially depen- 
dent lattice spacing to be the distance between the 
lattice sites k and k + 1. Our interaction coefficients are 
then given by 



Vk,k+l — 



E i= i a 



(k+j-i) 



(61) 



A. Alternating Spacings 

We initially consider the specific case of a lattice with 
two alternating spacings <2i, 02, such that the lattice sites 
are effectively grouped in pairs separated by the smaller 
spacing. For an even number of lattice sites this pre- 



serves the reflection symmetry of the Hamiltonian. For 
simplicity we introduce the ratio T = a\/a% between the 
lattice spacings. Then, the next-neighbour interactions 
are given by V := C 6 /af and W := C 6 /a% =T 6 V, re- 
spectively. For Vt/V < 1 we can again realize a weak 
laser coupling regime that allows for an analytical inves- 
tigation similar to Sec. |III A| We split up the Hamilto- 
nian H = Hq + Hl + .H/nt m t° a dominant contribution 
Hq = + Hd given by the familiar laser detuning op- 
erator and the next-neighbour Rydberg interactions and 
a perturbation given by the comparatively small laser 
coupling Hamiltonian and the long range Rydberg inter- 
actions Hj£ t . For an even number of lattice sites, the 
next-neighbour interactions may be written as 



N/2 



N/2 



n m n (2k+i) 



fc=i 



fc=i 



= v(nZ + tnZ) 



Here, we have already divided the next-neighbour inter- 
actions H™ t into pairs separated by a\ and 02, respec- 
tively, by introducing the corresponding next-neighbour 

NY. = y^ni 2fc -V 2fc) and N™ = 



pair operators 

v-»JV/2 ^(2fe)^(2fe+l) 



respectively. In the above-mentioned 
case of a lattice with an odd number of sites, these two op- 
erators must be appropriately modified, but most of the 
general spectral features are similar. The unperturbed 
Hamiltonian Hq is thus given by 



H = A (N e - N/2) + V [NY e + T 6 N™] 



(62) 



Since Hq is diagonal in the canonical product basis, this 
again leads to a simple formula for the eigenvalue of Hq 
for a canonical basis state | S) £ Su , 



E T (S) = A[N e (S) 



N/2]+V[Nl(S)+r e N^(S)], 

(63) 



where NY e (S) and N^ (S) are the eigenvalues of the ex- 
citation pair operators N^ e and for the state \S). 



Figure 10 presents the resulting energy spectrum and 
excitation pattern of the ground state obtained for T 6 = 2 
and an odd number of lattice sites N = 13. Comparing 
this result with the case of a constant spacing (equiva- 
lent to T = 1) shows a qualitatively different excitation 
pattern for the energetic ground state for A < —V. For- 
mally, the resulting energy eigenvalues [cf. ( |63[ )] share 
some similarities with our result obtained for a constant 
lattice spacing [cf. (fl2|]. V plays the role of V\ while 
:= [NY e (S) + r r N™{S)] may be understood as an 
effective, in our case of T 6 = 2 integer, pair number. It 
is clear that, depending on the specific value of T 6 , this 
generalization allows for more kinds of degeneracies as 
observed in figure 10 'a). Although there exist further 
combinatorial constraints on (S) and NY e (S) (e.g., 
^ee + -^ee = N ee ) it is now possible to bring two states 
featuring the same number of N e excitations but differ- 
ent excitation pair numbers into degeneracy. This cannot 
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FIG. 10: (a) Spectrum and (b) local excitation probabil- 
ity for the energetic ground state in a patterned lattice with 
N = 13, T = ^ = v^2 and fi = 0.05V. In contrast to the 
case of a single lattice spacing (cf. figure nj) we see that for 
A G [— 2V, — V] the ground state is given by superpositions 
of states featuring different excitation patterns, especially for 
A/V < — 1. In panel (b) white regions correspond to unity 
Rydberg excitation probability. 



be achieved dynamically, as this would require a varia- 
tion of the lattice constants during the experiment, which 
is likely to destroy the quantum coherence of the sys- 
tem. But it may be exploited to adiabatically prepare 
superpositions of lattice states with different pair num- 
bers N^,NY e by means of a time-dependent detuning 
A(t). 

We also remark at this point that the results obtained 
so far already have an interesting consequence: distinct 
spectral features will depend in some fashion on the sixth 
power of the ratio of the lattice spacings T 6 . This ex- 
ponent is due to the assumed van-der-Waals interaction 
potential which scales with the interatomic distance R as 
R~ 6 . For a general scaling V(R) cx R~ n the ratio would 
enter as F n . Hence, an experimental study of such a sys- 
tem could exploit this to verify the spatial dependence 
of the interaction potential. This would be achieved by 
choosing values for A, V and T which lead to very specific 




-2 2 
Site position x[ao] 

FIG. 11: The crystal states | [N e ]) (excitations marked in 
red) are presented for N e = 2, 3, 4, 5, 6, 7 in a lattice of size 
No = 61. As can be seen, to realize these states, only N = 13 
lattice sites (lowest row) must actually be occupied by atoms. 



and ideally experimentally testable degeneracies. How- 
ever, this goes beyond the scope of the present work and 
will be subject of a separate study. 



B. Crystal State Transitions 

A further interesting sequence of ground state transi- 
tions follows from our previous investigation of the suc- 
cession of crystal states in the full blockade regime. For 
N = 13 we can realize pronounced transitions |[2]) — ► 
| [3]) — > | [4]) — > | [5]) for the case of a single-spaced lattice 
(see section III I. For a lattice of size N = 61, on the 
other hand, we would be able to achieve two additional 
transitions | [2]> -> | [3]> -> |[4]) -> \{5}) -> | [6]) -> \[7}). 



-> 

However, the energy spacing between the crystal states 
scales with (N — 1)~ 6 for a fixed number of excitations, 
rendering an experimental realization impractical. This 
drawback can be compensated, though, by decreasing the 
lattice constant an. 

There exists an alternative possibility of realizing the 
N = 61 crystal state transitions with a strongly non- 
uniform lattice consisting of only N = 13 lattice sites. 
In principle this can be thought of a full N = 61 lat- 
tice, where the atoms have been removed from 48 sites. 
Figure[TT]illustrates the corresponding spatial pattern ac- 
cording to which the sites have to be arranged. For the 
case of figure [TT] the values of the detunings at the tran- 
sition points are given by (25) with N — 13, providing 



a more favourable scaling of the energy spacing between 
the crystal states. As indicated above, the same tran- 
sition detunings can be achieved by a uniform TV = 61 
lattice with a reduced lattice spacing of 0.2 an. Here, 
an denotes the spacing of a uniform lattice consisting of 
N = 13 sites. The redistributed N = 13 lattice, on the 
other hand, only demands a minimal spacing of 0.4 an. 
If the minimal lattice spacing is experimentally limited, 
the redistributed N = 13 lattice thus provides a more 
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FIG. 12: (a) Spectrum and (b) local excitation probability 
for the energetic ground state in the patterned lattice of figure 
11 with fl = W~ 4 V. As discussed in the text, this particular 
setup emulates a iVo = 61 site lattice. The transitions between 
the crystal states | [N e ]) are evident in panel (b) where white 
regions correspond to unity Rydberg excitation probability. 



weak laser regime, a systematic classification of the state 
space has been performed. In the so-called full block- 
ade regime, where no neighbouring Rydberg excitations 
in the lattice are allowed, the transitions between crystal 
Rydberg lattice states were identified. Turning our inves- 
tigations to the strong laser regime, an extensive analytic 
treatment allowed the full characterization of the spec- 
trum by performing site-specific rotations in spin space 
and introducing fermionic ladder operators. An excellent 
agreement of our perturbative treatment and the numeri- 
cally obtained spectrum has been found. Going from reg- 
ularly spaced to patterned lattices, we investigated two 
particular examples. For alternating lattice constants, 
an additional energy scale is introduced that alters the 
spectral properties as well as the local Rydberg excita- 
tion probability of the ground state. In addition, we have 
highlighted a specific setup involving an irregular lattice 
site distribution that allows the emulation of the crys- 
tal state transitions of a more extended one-dimensional 
lattice. 
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Appendix A: Perturbative Energy Corrections in 
the Strong Laser Regime 



favourable scaling of the relevant detunings. 



Figure 12 a) presents the energy spectrum with the 



first five transitions between crystal states. In panel (b) 
the corresponding local Rydberg excitation probability 
is illustrated. One can nicely observe the transitions be- 
tween the crystal states. In order to get the sharp tran- 
sitions between the crystal states, a relatively small Rabi 
frequency of £1 — 10~ 4 V is assumed in this figure. Due 
to the sensitive scaling of the corresponding transition 
detuning values with the interaction exponent n, their 
exact experimental measurement could provide a useful 
test of the validity of the assumed interaction potential. 



Here we present the formulas used for calculating 
the perturbative energy corrections presented in Section 
|III B| This is achieved by re-expressing the perturba- 
tion operators in terms of the fermionic ladder operators 

Vk = YliLiRklCi- The coefficient matrix (Rkl)k,i=i is 
real, symmetric and orthogonal: 



Rki — Rik = 

N N 



2 s i n kl7T 
N+l Bm N+l ' 



^ RklRjl = ^ RklRlj — $kj- 



1=1 



BRIEF SUMMARY 



This can be proved by means of the following identities. 
We find that for m G Z 



Focusing on the frozen gas regime at ultracold temper- 
atures, we have provided a comprehensive analysis of the 
spectral properties of the coherent, laser-driven Rydberg 
excitation in one-dimensional ordered systems. In the 



N 



^2 exp(i 



It Itl K ' 

N+l' 



71=0 




sin( 



JV+1 ' 



for m G 2(N + 1)Z, 



l+COs( ivT-Ll ) 

1 + i , mY for odd m, 



otherwise. 
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Prom this, it is straightforward to show that for k, I € 
{1,2,. ..,N} 

N 

1 

~ ~ 2 



N 



exp|z 



n=0 



n{k—l)n i 1 

W+i if 



N+l 



>ki- 



and one finds therefore 

N 

E sin (^)cos(^) 



n=l 



1 N 

2 3 H { exp [ ? 



n(fc+/)7r I 

JV+T J 



+ exp[i 



n=0 



i(k— 1 
V+l J J 



. , k7T . 

sm (]y+T) 
cos (]V+T) _cos (lV+T' 



JV+1 

for even (k ± Z) 
for odd (k ± 



Where 5ft, 5 indicate the real and imaginary part. The 
unperturbed XY-Hamiltonian reads 



H 



xy 



JV 

= X] VkVk 

k=l 



n 



Vi , kir , 

T cos( ]vTi' 



const. 



(Al) 



and the perturbations are given by the boundary term 
Vi 



H h = 



r(l) 



and the off-resonant couplings 



H' 2 = 



v N ~ 1 

T ^ 

fc=i 



ai fe) ai fe+1) 



(A2) 



(A3) 



The first energy manifold is spanned by the ground state 
G\. The first excited energy manifold is given by the 
set of states that contain a single excitation 



G 



1< k < N) 



One finds that 



H 



■j [4 + ci + (c 



1/ N 

fc=i 



-Rife 



1 + ( _ 1) fe+JV + 



l-(-l) 



Since the last term in ( A5 1 is antisymmetric with respect 
to k and Z and they are summed over their full range, it 
is convenient to calculate the antisymmetrized value of 
the bracketed expression. We define thus 



■ kl ■- 



'N—l 



N 

N+l X 

n=0 

N 



Rk,nRl,n+lj — (k *H> I) 

{sin(^)si 



r (n+l)/7ri 
L N+l ' 



(k++l)} 



NTi E { H#TT) cos( 1 ^ T ) 

n=0 

+ s i n (^L)cos(^ I )sin( 1 ^ T )] - (k ++ I)} 
for even (k ±1), 

for odd (k±l). 



■'(1 



-)sin(- 



cos (jv+T 



, fcir , 

- c ° s (jv+t) 



The off-resonant couplings are then given by 



V N 

k,l = l 



We thus see explicitly that the perturbation operators 
only couple manifolds of different numbers of fermionic 
excitations, specifically AJV + = ±1 for Hf, and AiV + = 
±2 for H'2, respectively. This implies that all first or- 
der energy corrections vanish due to the fact that N + is 
diagonal in our basis. In our specific case we find that 
we can calculate the second order contributions due to 
Hf, and H' 2 independently, since for any given eigenstate 
these operators couple to different subspaces. The rele- 
vant matrix elements for the corrections to the ground 
state are given by 



G 



VkH b 



G) = ^R lk [l-(-l) k ], 



G 



mviH'2 



G 



Vx 



Ft 



kl- 



From these results and the specific form of Fki we see that 
only specific states contribute to the energy corrections. 
This is due to the overall symmetry under reversal of the 
lattice, which is preserved by the perturbations. For the 
Vk J corrections to the states of the first excited manifold we 
need the following matrix elements 



where we have used R.Nk — {~l) k+1 Rik- The second 
part of the perturbation is 



H'o = 




n c n+l c n c n+l 



Rk,nRl,nA 



VkVl 



(A4) 
(A5) 



G 



G 



r]kmHb% g 



VkVlVmH^vl 



G 



(FklS 7nq — Fk m Sl q + Fi m Skq) ■ 



1G 



The ground state energy correction is thus given by Making the same additional approximation as above, we 



G G,b ^ G,off ' 



find for the corrections to an excited state 77+ 



2 



AE 



T/2 E>2 
(2) \ - -"-1 



G.i, iq X! 1 , ; , . .... ■ .:/.■ un 
fc=i 1 + 2S COb JV+i 



,2fe-l 



x/2 N N 



fc=l Z=fc+1 



G 



A £( 2 ) — ^ f f ^ A^») = AE^ + AE% + AE^ 

Ai W*^.* 2 ' 

can further approximate the energy corrections by ne- q '° 4fj 

glecting the next to leading order contributions in the 
denominators: 



Before we continue by writing down the energy correc- 
tions to the singly excited states we should note that we AE^ ~ ^ R 2 ^ ^~ ^* 
can further a.nnroximate the enerev corrections bv ne- 1> b ACi ^-^ 9 



fc=l 



V 2 N / JV 



2fe-l) 



vl 
8n 
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